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Abstract: This five-year study (2016-2021) in Morocco’s Mediterranean climate investigated the 
effect of nitrogen fertilization and genotypic selection on wheat yield and associated components. 
Utilizing a split-plot design, the study assessed five wheat genotypes— Faraj’, ‘Luiza’, ‘Itri’, ‘Karim’ 
and ‘Nassira’—under three nitrogen application rates (120, 60 and 0 kg/ha) across thirty plots 
with two replicates. Interactions between nitrogen and year showed marked significance in yield 
(p = 0.001), biomass (p = 0.002), TKW (p = 0.003) and Spk/ m? (p = 0.001), underscoring the variability 
in optimal nitrogen application rates across different years. Additionally, significant interactions 
between variety and year were observed for biomass (p = 0.001) and G/ m? (p = 0.001), indicating 
variability in the performance of different varieties across years. The ‘Itri’ genotype showed the 
highest yield in 2017, while ‘Luiza’ was pre-eminent in 2018, with ‘Ttri’ producing the most biomass. 
‘Faraj’ demonstrated consistent superiority in yield and biomass during 2019 and 2020. Our integrated 
principal component analysis and quadratic models elucidated that an intermediate nitrogen rate of 
60 kg/ha (N2) was particularly advantageous for the ‘Faraj’ and ‘Karim’ genotypes. These findings 
highlight the substantial impact of informed nitrogen level adjustment and genotypic selection on 
yield optimization. 
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1. Introduction 


Wheat and barley, the predominant cereal crops in the Mediterranean region, are 
paramount for regional food security [1]. In the context of the Mediterranean region, 
optimized cultivation of wheat has far-reaching implications not only for food security but 
also for sustainable agricultural practices [2,3]. 

In Morocco, wheat accounts for 65% of the country’s total cereal area, covering approx- 
imately 2.9 million hectares and making it the most significant crop [4,5]. Durum wheat 
(Triticum turgidum var. durum Desf), an ancient wheat species, forms a substantial part of 
human and animal diets [6]. However, despite the crop’s economic importance, durum 
wheat grain yields in Morocco remain low and unpredictable in dry areas due to limitations 
such as low and poorly distributed rainfall [7-10]. 

Research in the field has highlighted challenges in wheat yield improvement, notably 
the inverse relationship between cereal protein content and grain yields, which compli- 
cates genotypic selection strategies for breeders [11]. This negative correlation implies 
that selecting higher yields often results in lower protein levels in the grains, presenting 
a complex balance in genotypic selection strategies [12]. Yield at harvest is a genetically 
controlled trait [13,14]. However, environmental factors and farming practices often lead to 
variations in yield [15-17]. Yield variability in wheat is influenced by several agronomic 
factors, with nitrogen management being crucial among them. Proper nitrogen fertilization 


Nitrogen 2024, 5, 266-287. https: / /doi.org/10.3390/nitrogen5020018 


https: //www.mdpi.com/journal/nitrogen 


Nitrogen 2024, 5 


267 


can significantly enhance yield and biomass production by optimizing the plant’s nutri- 
tional status, especially in regions challenged by environmental constraints and resource 
limitations [10]. This approach aligns with our study’s focus on assessing the impacts of 
nitrogen fertilization and genotype selection on wheat productivity. 

To ensure optimal yield, it is crucial to identify farming practices, which cater to 
this objective. One such practice is the strategic management of nitrogen fertilization, 
tailored to the crop’s needs at various developmental stages [18-20]. By understanding the 
relationship between nitrogen absorption and grain yield [21], farmers can better manage 
nitrogen inputs according to the specific requirements of the chosen variety [22]. 

With the increasing global demand for wheat, developing varieties with improved 
productivity and resource utilization is vital. Therefore, this study aims to assess the effects 
of nitrogen fertilization and genotype selection on wheat yield, biomass production and 
yield components. It seeks a comprehensive understanding of the intricate interdependen- 
cies among these variables. Through a five-year study (2016-2021), we specifically focused 
on comparing the performance across three categories of wheat varieties: new (‘Itri’ and 
‘Luiza’), medium-aged (‘Faraj’ and ‘Nassira’) and ancient (‘Karim’). This comparison aimed 
to examine the influence of nitrogen application rates on each category to identify optimal 
strategies for maximizing wheat productivity. 


2. Materials and Methods 
2.1. Experimental Sites and Climate Conditions 


This study was conducted over five years, from 2016 to 2021, at three experimental 
stations belonging to Morocco’s National Institute of Agricultural Research (INRA). 

Table 1 provides detailed geographic coordinates, ecosystem types and five-year 
precipitation trends of these INRA experimental stations in our study, encapsulating key 
environmental data from the period 2016-2021. 


Table 1. Geographic coordinates, ecosystem types and five-year precipitation trends of INRA experi- 
mental stations in our study (2016-2021). 


. . 2016/2017 2017/2018 2018/2019 2019/2020 2020/2021 
Station Location Ecosystem bsg gs es igsogs Sevag Renee a ds 
Precipitation Precipitation Precipitation Precipitation Precipitation 
Name (Lat, Long) Type 
(mm) (mm) (mm) (mm) (mm) 
SidiEl Aydi 33.12218° N, Semi-Arid, 
(SEA) 7.63315° W Rainfed ae ae ae Ber sad 
Merchouch 33.61375° N, — Favorable/Sub- 
(MCH) 6.71786° W Humid = saa ae as a 
Tassaout 31.82021° N, | Dry and Hot, 
(TST) 7.43806° W Irrigated saa aD — ad ce 


The locations of these stations and regional variations are depicted in Figure 1, which 
presents dual maps illustrating the growth season duration (in days) and annual precipita- 
tion (in mm/year) across various study locations and regions in Morocco for the 2016-2017 
to 2020-2021 agricultural seasons. 
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Figure 1. Dual maps of growth season duration (in days) and annual precipitation (in mm/year) 
across different study locations and regions in Morocco, spanning the 2016-2017 to 2020-2021 


agricultural seasons. 


2.2. Experimental Design 

The experiment was conducted using a split-plot design with two factors: nitrogen 
treatment (N) in the main plots and variety (V) in the sub-plots. Three levels of nitrogen 
(N) factor were used: N1 = 120 kg/ha, N2 = 60 kg/ha and N3 = 0 kg/ha. Five different 
varieties were used in the sub-plots. The experiment was repeated twice, resulting in a total 
of 30 plots. The characteristics of the varieties are summarized in Table 2. 


Table 2. Variety, year of registration, drought stress response and Hessian fly reaction. 


Variety Year of Response to Stress: Hessian Fly Pedigree 
Release Drought Resistance 
‘Faraj’ 2007 Moderate resistance R Hybrid Nassira, Qarmal, Lahn (ICARDA) 
Ter’ 2017 Tlenaat S RISSA/GAN//POHO_1/3/ PLATA_3 / /CREX/ALLA/x 
Karim 

‘Karim’ 1985 Moderately tolerant S Bittern ‘S’ « JO’S’. AA”:S’/ /FG’S’ » 

‘Luiza’ 2011 Moderate resistance S RASCON_39/TILO_1 
; es INRA Selection on CIMMYT EI, 12 TA14/BD3/ /Isly # 
Nassira 2003 Tolerant R CF41530-1548 
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2.3. Plot Characteristics and Agricultural Practices 


Each experimental plot occupied a 2.7 m2 area, with dimensions of 2.5 m in length 
and 1.08 m in width. To ensure optimal growing conditions, standard agronomic practices, 
such as soil preparation and weeding, were employed throughout the growing season. 
Ammonium nitrate fertilizer with a 33.5% nitrogen content was applied in two stages 
to align with the crop’s growth phases: half at the tillering stage and the remaining half 
at the stem extension stage, ensuring targeted nutrient support. Seeds were sown at 
a density of 400 g/m? around mid-November using a Wintersteiger plot seeder with an 
inter-row spacing of 20 cm and a sowing depth of 3.5 cm to accommodate the planting 
dates recommended for the three different regions under study. To manage pests and 
avoid the peak cycle of the Hessian fly, seeding dates were carefully selected, and manual 
weeding was conducted to control weed competition. Harvesting was completed by the 
end of June for all study years, marking the conclusion of the annual growth cycle. 


2.4. Data Recording and Sampling 


Data were collected from three distinct locations on six key agronomic traits, adhering 
to rigorous methodologies and international standards. Yield was initially measured in 
¢/2.7 m? and subsequently converted to kg/ha for standardization. Biomass, quantified 
pre-harvest, included only the above-ground plant material and was weighed in kilograms. 
The thousand-kernel weight (TKW) was determined using an electronic grain counter com- 
pliant with NF V03-702 [23] and ISO 520 [24] standards. The number of spikes per square 
meter (Spk/m/7) was assessed based on designated one-square-meter sample areas within 
each plot. 


2.5. Data Processing 


Data analysis involved multiple software tools. In Excel, we computed the mean of the 
data collected from the three locations for later calculations. For graphical representations, 
such as bar plots and GT bi-plots, we used R. 

Minitab 18 was specifically employed for statistical analysis. A three-way ANOVA 
was conducted to understand the effects of variety, nitrogen, year and their interactions on 
the observed measurements. 

Lastly, OriginPro was utilized for principal component analysis (PCA). PCA aided in 
simplifying data visualization and interpretation by reducing the dimensionality of the 
dataset while retaining variance information. 


3. Results 
3.1. Genotypic Influence on Yield and Yield Components 


Initial analysis of the ANOVA results for yield and biomass demonstrated significant 
variation among the durum wheat varieties studied. Both factors displayed a significance 
level of p < 0.05. This finding emphasizes that the observed differences in yield and biomass 
among the varieties are not due to chance alone but likely reflect distinct characteristics 
inherent to their genotypes. 

During the five-year research span (2016-2021), a trend was observed where newer 
varieties tended to achieve higher yields than the oldest variety. Among the evaluated 
varieties, ‘Faraj’ exhibited superior productivity in yield (3430 kg/ha) and biomass produc- 
tion (12,196 kg/ha). Meanwhile, ‘Luiza’ and ‘tri’ also exhibited notable yield and biomass 
production performance. ‘Luiza’ demonstrated a yield of 3422 kg/ha and a biomass of 
11,949 kg/ha, while ‘Itri’ produced a yield of 3347 kg/ha and a biomass of 11,889 kg/ha. 
These results indicate the potential of both new and moderately new varieties in agri- 
cultural production. The older varieties, such as ‘Karim’ (yield: 3239 kg/ha, biomass: 
11,267 kg/ha) and ‘Nassira’ (yield: 3261 kg/ha, biomass: 10,543 kg/ha), demonstrated 
lower but competitive productivity. The average yield and biomass values calculated across 
all varieties stood at 3340 kg/ha and 11,569 kg/ha. However, it is important to note that 


Nitrogen 2024, 5 


270 


while the ANOVA results indicated significant yield differences among the varieties, the 
effect of variety on biomass was not statistically significant (ns). 

The study also explored grain size using the metric of TKW. The ANOVA analysis 
for TKW revealed a highly significant difference among the varieties, denoted as p < 0.001. 
Regarding TKW, ‘Karim’ emerged as the leader with the greatest value of 42.46 g. ‘Ttri’, 
however, recorded the lowest TKW of 39.21 g. 

A similar trend was observed with Spk/ m2? (plant number per square meter). The 
ANOVA analysis indicated significant differences among the varieties at p < 0.01. Here, 
‘Faraj’ once again asserted its superiority, recording the greatest number of spikes per square 
meter (332 Spk/m/7). In contrast, ‘Karim’ had the lowest Spk/m2, registering 294 Spk/m?. 
The average plant density across all the varieties was 311 spikes per square meter. 

Regarding G/ m? (grain number per square meter), the ‘Luiza’ variety proved to be the 
most prolific, producing the highest number of grains per square meter at 7784. ‘Nassira’, 
on the other hand, had the lowest grain count at 7152 grains per square meter. Across 
all varieties and treatments, the overall mean G/m2 was determined to be 7460 grains per 
square meter. 

The research also evaluated the characteristics of grains per spike (G/S), showing 
significant differences among the varieties at p < 0.05. ‘Faraj’, ‘Luiza’ and ‘Itri’ were the 
top performers, each producing 27 grains per spike. Conversely, ‘Karim’ and ‘Nassira’ 
registered the smallest G/S values, with 26 and 24 grains per spike, respectively. These 
findings highlight the distinct variation in G/S across the different durum wheat varieties. 
The variations observed in this study can be attributed to the distinct genotypes of the 
varieties tested, suggesting that genetic factors play a pivotal role in determining these 
agronomic traits. 

A visual representation in the bi-plot analysis (Figure 2) clearly presents the relation- 
ship between the wheat varieties and agronomic traits. The positioning of each variety 
on the bi-plot provides insights into their respective agronomic characteristics. The new 
varieties, ‘Luiza’ (2011) and ‘Itri’ (2017), showed a close association with grains per square 
meter (G/m?) and the number of grains per spike (G/S). Conversely, the ‘Nassira’ (2003) 
variety was associated more closely with the thousand-kernel weight (TKW). The ‘Faraj’ 
(2007) variety was more closely aligned with spikes per square meter (Spk/m7?) than with 
biomass or yield. Meanwhile, the older variety ‘Karim’ (1985) found itself on the opposite 
side of the bi-plot, indicating differing agronomic characteristics. 


GT Biplot 
Scaling = 1, Centering = 2, SVP =2 


Spk/m2 
Faraj 
Nassira 


Biomass (q/ha) 


PC2 (31.1%) 
g 


Karim 


i) 
PC1 (58.52%) 


Figure 2. Genotype by trait (GT) bi-plot. The trait codes are as follows: TKW = Thousand-kernel 
weight, Spk/ m2 = Spikes per square meter, G/ m2 = Grains per square meter, G/S = Grains per spike. 
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3.2. Influence of Nitrogen on Biomass, Yields and Yield Components 


Turning to the role of nitrogen, the ANOVA results showed a significant effect on 
crop yield (p < 0.05), biomass (p < 0.001) and thousand-kernel weight (TKW) (p < 0.001), 
as outlined in Table 3, demonstrating nitrogen’s significant influence on these metrics. 
The maximum yield was obtained following an application of 120 kg/ha of nitrogen, 
reaching 3404 kg/ha, while the minimum yield was observed with no nitrogen application, 
at 3234 kg/ha. 


Table 3. Significance levels of ANOVA for genotype and nitrogen effects on durum wheat yields, 
yield component. 


Source of Variation eyika) “igh TKW (g) Spk/m? G/m? G/S 
Variety 

‘Faraj’ (2007) 3430 * 12,196 # 41.28 > 332 * 7484 * 27? 

‘Itri’ (2017) 3347 4 11,889 a> 39.21 © 305 > 7422 4 274 

‘Karim’ (1985) 3239 > 11,267 > 42.46 * 294. 7457 * 26° 

‘Luiza’ (2011) 3422 2 11,949 @ 39.98 © 302 > 7784 * 27% 

‘Nassira’ (2003) 3261 10,543 ab 42.32 4 311 2 71524 244 

Mean 3340 11,569 41.05 311 7460 264 

Nitrogen 
N1 120 kg/ha 3404 4 12,266 @ 41.68 * 305 ab 7352 * 25 4 
N2 60 kg/ha 3376 * 12,2774 40.56 » 319 4 7730 * 27° 
& 

N3 0 kg/ha 3234 > 10,917 > 40.92 > 302° 7549 @ 274 

Mean 3338 1182 41.05 308 7544 26 

ANOVA DE eyika) ‘ugha = TKW(@) Spin? G/m? G/s 
Variety 4 * ns AYE is ns . 
Nitrogen 2 id a wee * ns ns 
Year 4 EE HAE HE Cee era ated 
Variety x Nitrogen 8 ns ns ns ns ns ns 
Variety x Year 16 EH AK EE See ra ns 
Nitrogen x Year 8 — = ~ sil * ns 
Variety x Nitrogen x 32 ng ie ne He ns ce 


Year 


TKW = Thousand-kernel weight, Spk/m? = Spikes per square meter, G/m? = Grains per square 
meter, G/S = Grains per spike. . ANOVA—Analysis of variance, where * = Significant at p < 0.05; 
* = Significant at p < 0.01; *** = Significant at p < 0.001; ns = Not significant. Means without a common let- 
ter also vary greatly. 


Regarding biomass, the ANOVA analysis also revealed a significant impact of nitrogen 
(p < 0.001). Consistent with the yield results, the greatest biomass was achieved with 120 kg 
of nitrogen, reaching 12,277 kg/ha. In contrast, the smallest biomass was noted for 0 kg/ha 
of nitrogen, standing at 10,917 kg/ha. 

Additionally, the ANOVA analysis detected a significant effect of nitrogen on the TKW 
of the crop (p < 0.001). In alignment with this, the application of 120 kg/ha of nitrogen 
yielded the greatest TKW, at 41.68 g. The smallest TKW, however, was noted with 60 kg/ha 
of nitrogen, registering 40.56 g. 

The analysis of variance (ANOVA) findings indicated a significant difference in the 
impact of nitrogen on the quantity of spikes per square meter (Spk/m?) with a p = 0.003. 
Notably, utilization of 60 units of nitrogen had the most pronounced effect on the spike 
count per square meter, resulting in a total of 319 Spk/m?*. On the other hand, the ab- 
sence of nitrogen application, marked as 0 units, yielded the lowest spike count, totaling 
302 spikes per square meter. 

Furthermore, the application of nitrogen showed no significant impact on the quantity 
of grains per square meter (p = 0.095) or the number of grains per spike (p = 0.193). The 
experimental results indicated that the maximum grain yield per unit area resulted from 
the precise application of 60 units of nitrogen, which resulted in a grain count of 7730 per 
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square meter. Conversely, the minimum grains per square meter were observed with the 
application of 120 kg/ha of nitrogen, amounting to 7352 grains per square meter. The mean 
grains per square meter stood at 7544. Similarly, the highest number of grains per spike 
occurred with 0 kg/ha and 60 kg/ha of nitrogen, both delivering 27 grains per spike. The 
lowest number of grains per spike was obtained with 120 kg/ha of nitrogen, registering 
25 grains per spike. 

The selection of nitrogen doses of 120 kg/ha and 60 kg/ha was informed by a com- 
prehensive review of preliminary trials and relevant literature, aiming to identify optimal 
rates, which balance environmental sustainability with yield improvements. 

A closer examination of nitrogen use efficiency (NUE) reveals a more nuanced picture 
of nitrogen’s role in wheat production. Specifically, the NUE for the 60 kg/ha dose stood 
at 2.37 kg grain/kg N, surpassing the 1.42 kg grain/kg N efficiency observed for the 
120 kg/ha dose. This differential efficiency highlights the diminishing returns on yield 
with increased nitrogen application and underscores the importance of optimizing nitrogen 
management to achieve sustainable wheat production. 

The findings from this study illuminate the critical balance between nitrogen appli- 
cation and wheat yield optimization. While higher nitrogen doses can enhance yield and 
biomass, the efficiency of nitrogen use, as evidenced by the NUE calculations, suggests 
that lower application rates may offer a more sustainable pathway for maximizing wheat 
productivity. This balance is crucial for agronomic practices, especially in regions where 
environmental sustainability and resource optimization are paramount. 

As illustrated in Figure 3, nitrogen application rates and variety significantly influence 
crop yield. Recent varieties (‘Faraj’, ‘Itri’ and ‘Luiza’) exhibit greater yields compared 
to older varieties (‘Karim’ and ‘Nassira’) at a given nitrogen application, which can be 
attributed to advancements in plant breeding. Furthermore, nitrogen application is pivotal 
in yield enhancement, with N1 (120 kg of nitrogen) and N2 (60 kg of nitrogen) consistently 
leading to superior yields over N3 (0 kg of nitrogen) across each variety. The variability in 
response to nitrogen levels among the varieties underscores the differences in nitrogen use 
efficiency and yield potential inherent to each genotype. 


Mean Yield for each Genotype-Nitrogen combination 


3572 
2451 si 3504 2442 
3330 wea 
19 
2268 3224 3245 2288 4104 3248 570 

3000 
Nitrogen 

2000 
3] N1 120 kg/ha 
| N2 60 kg/ha 
BB nz okona 

1000 

0 


Faraj ITRI Karim Luiza Nassira 
Genotype 


Mean Yield (kg/ha) 


Figure 3. Yield response to nitrogen doses across different varieties. 
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3.3. Insights on Interaction Effects of Agronomic Traits 


Expanding the scope of the analysis, the study explored the interaction effects among 
variety and nitrogen, variety and year, and nitrogen and year, as detailed in Table 3. 

Notably, the relationship between variety and nitrogen showed no significant dif- 
ferences in several agronomic traits—yield (p = 0.909), biomass (p = 0.475) and Spk/m? 
(p = 0.954)—according to Table 3. These data imply that the influence of nitrogen re- 
mains stable across various varieties. This means that the outcome of nitrogen application 
is likely to be uniform irrespective of the wheat variety chosen, streamlining nitrogen 
management decisions. 

When examining the interaction between variety and year, significant effects were 
observed for yield (p < 0.001), biomass (p < 0.001), TKW (p < 0.001), Spk/ m2 (p < 0.001) 
and G/m? (p < 0.001), but not for G/S (ns), indicating that the performance of different 
wheat varieties can vary significantly from year to year, depending on environmental 
conditions and other factors. Similarly, the interaction between nitrogen and year showed 
significant differences for yield (p < 0.001), biomass (p = 0.002), TKW (p = 0.003) and Spk/m? 
(p < 0.001), with a noticeable effect on G/m? (p < 0.05), highlighting how optimal nitrogen 
application rates may vary across different years, influenced by annual weather patterns 
and other environmental factors. Such variations are pivotal to consider in agronomic 
planning. The three-way interaction between variety, nitrogen and year was predominantly 
found to be non-significant, pointing toward a stable combined effect of these parameters 
on the agronomic traits. 

Venturing into a deeper data analysis, the principal component analysis (PCA) (Figure 4) 
unravels the associations between traits and nitrogen application rates. The principal 
component analysis (PCA) of the data reveals two principal components, PC1 and PC2, 
which account for 57.71% and 42.29% of total variance, respectively. The TKW is closely 
associated with the greatest nitrogen application rate (N1, 120 kg/ha), indicating a positive 
correlation between these two variables. Similarly, the number of spikes per square meter 
and grains per square meter (G/m7) are closely associated with the N2 (60 kg/ha) nitrogen 
application rate. In contrast, yield and biomass are situated on the opposite side of the PCA 
plot, indicating a negative correlation with the N3 0 kg/ha application rate. Interestingly, 
the yield and biomass are situated between the N1 and N2 nitrogen application rates on 
the PCA plot. This suggests that an intermediate nitrogen application rate, between N1 
and N2, could potentially result in optimal yield and biomass. 


-1.0 


PC 2 (42.29%) 


=i5 -1.0 -0.5 0.0 05 1.0 
PC 1 (57.71%) 


Figure 4. Principal component analysis (PCA). The trait codes are as follows: TKW = Thousand-kernel 
weight, Spk/ m2 = Spikes per square meter, G/ m2 = Grains per square meter, G/S = Grains per spike. 
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To further explore the relationship between nitrogen application rates and yield, we 
employed quadratic models for two distinct crop varieties, ‘Faraj’ and ‘Karim’. These 
models were derived from extensive field trials and were validated to ensure an accurate 
fit with empirical data. The equations are as follows: 

e For the ‘Faraj’ variety (Figure 5): Yparaj(x) = —0.0596x? + 8.690x + 3265.4. 
e —- For the ‘Karim’ variety (Figure 6): Yxarim(X) = —0.0204x? + 2.97x + 3183.4. 


Yield vs. Nitrogen Rates for Faraj Variety 


3600 - Y=- 0.0596x? + 8.690x + 3265.4 


3500 - 


ic 
€ 

® 
=3400- 
Bf 

2 
2 

3300 - 
N 
3200 - 
0 0 100 
Nitrogen Rates (N) 
Figure 5. Quadratic response of yield to varying nitrogen rates in ‘Faraj’ variety. 
Yield vs. Nitrogen Rates for Karim Variety 
S200 = Y =-0.0596x? + 8.690x + 3265.4 
om Nt 

g 

S 

=3200- 

2 

= N3 


3100 - 


0 50 100 
Nitrogen Rates (N) 


Figure 6. Quadratic response of yield to varying nitrogen rates in ‘Karim’ variety. 


Where Y represents the yield in kg/ha, and x represents the nitrogen application rates 
in kg/ha. 
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Upon optimization, for the ‘Faraj’ variety, the nitrogen application rate corresponding 
to maximum yield was found to be x & 72.42 kg/ha. This nitrogen application rate predicts 
a maximum yield of approximately YFaraj © 3580 kg/ha. The second derivative confirmed 
this as a local maximum. 

Similarly, for the ‘Karim’ variety, the nitrogen dose for optimal yield was determined 
to be x & 74.25 kg/ha, with a corresponding yield of Yxarim © 3294 kg/ha. Again, the 
second derivative confirmed this as a local maximum. Although the optimal nitrogen doses 
for ‘Faraj’ and ‘Karim’ are nearly identical, indicating a similar nitrogen use efficiency, 
‘Faraj’ can achieve a higher yield than ‘Karim’ at this nearly equivalent nitrogen application 
rate, highlighting the potential for genetic factors to influence yield optimization even 
under similar nutrient management practices. 


3.4. Comparative Analysis of Wheat Varieties under Different Nitrogen Application Rates: 
A Five-Year Study on Yield, Biomass and Grain Characteristics 

In the inaugural year of our investigation (2017), we focused on assessing the impact 
of varying nitrogen application rates (N1, N2, N3) ona selection of wheat varieties (‘Faraj’, 
‘Ttri’, ‘Karim’, ‘Luiza’, ‘Nassira’), examining key parameters, including yield, biomass, 
TKW, number of spikes per square meter (Spk/m7?), grains per square meter (G/m7*) 
and the number of grains per spike (G/S). In this comprehensive study, a number of 
pivotal trends were identified. Across all nitrogen application rates, the ‘Itri’ variety 
demonstrated superior performance by achieving the greatest average yield of 3993 kg/ha, 
which corresponded to a notable 6.2% increase compared to the overall mean yield of 
3760 kg/ha, as illustrated in Table 4. In contrast, the ‘Karim’ variety exhibited the smallest 
average yield of 3502 kg per ha (kg/ha), representing a reduction of approximately 6.9% 
compared to the mean yield across all varieties. In terms of biomass, ‘Faraj’ notably 
excelled at the N2 level, accumulating 14,900 kg/ha—a 11.6% increase over the general 
mean biomass of 13,347 kg/ha. ‘Karim’, on the other hand, recorded the smallest biomass 
of 11,863 kg/ha at the N2 level—a 11.1% decline from the overall mean. Interestingly, 
although ‘Karim’ lagged in yield and biomass, it exhibited the greatest TKW of 49.06 g at 
the N1 level—a 14.8% increase over the general mean TKW of 42.74 g. Conversely, ‘Luiza’ 
registered the lowest TKW of 39.08 g at the N2 level, which was approximately 8.6% below 
the mean, according to Table 4. 

In the 2018 continuation of our study, we observed marked differences in wheat 
variety responses to nitrogen levels (N1, N2 and N3). ‘Luiza’ consistently outshone its 
counterparts in yield, averaging 4794 kg/ha, indicating an 8% rise from the revised mean 
yield of 4437 kg/ha for that year. However, at N3, ‘Luiza’ yielded slightly less than ‘Karim’, 
as highlighted in Table 5. Contrarily, ‘Faraj’ and ‘Nassira’ displayed suboptimal results, 
with yields of 4199 kg/ha and 4133 kg/ha, which were 5.4% and 6.9% below the mean, 
respectively. When it came to biomass, both ‘Ttri’ and ‘Luiza’ were prominent, especially at 
the N1 level, amassing an average biomass of 14,001 kg/ha, marking a 20.4% increase over 
the average biomass of 11,624 kg/ha. However, notably, ‘Faraj’ yielded less at the N3 level, 
posting a biomass figure of 10,128 kg/ha—12.9% less than the mean. Interestingly, while 
‘Faraj’ achieved the highest TKW of 48.35 g at the N3 level, surpassing the mean TKW of 
45.80 ¢ by 5.6%, ‘Itri’ fell short in this measure, with the smallest TKW of 43.46 g at the N2 
level—5.1% below the mean—as reported in Table 5. 

In 2019, the ‘Faraj’ variety consistently produced the largest yields across all nitrogen 
levels, with an average yield of 3425 kg/ha, representing a significant 16.5% increase over 
the overall mean yield of 2939 kg/ha, as illustrated in Table 6. In contrast, ‘Nassira’ and 
‘Itri’ generally underperformed in yield, registering average yields of 2628 kg/ha and 
2733 kg/ha, falling 10.6% and 7.0% below the mean, respectively. Regarding biomass, ‘Itri’ 
produced the highest average at the N2 level, with 13,406 kg/ha, which was a notable 28% 
increase over the overall mean biomass of 10,474 kg/ha. In contrast, ‘Nassira’ reported the 
smallest biomass of 7820 kg/ha at the N3 level, trailing the mean by 25.3%. Interestingly, 
despite its smaller yield, ‘Nassira’ led in TKW at the N1 level, with 39.7 g, exceeding the 


Nitrogen 2024, 5 


276 


Nitrogen 


overall mean TKW of 37.8 g by a modest 5%. ‘Itri’, however, produced a small TKW at both 
the N2 and N3 levels, recording the lowest average of 34.59 g—a decline of 8.5% from the 
mean—as reported in Table 6. 


Table 4. Comparative study of wheat varieties under various nitrogen levels in 2017: yield, biomass 
and grain characteristics with ANOVA insights. 


Yield Biomass 


Year Rate Variety (kg/ha) (kg/ha) TKW (g) Spk/m? G/m? G/S 
‘Faraj’ 3911 abc 14,136 abed = 42.97 bedef 4754 9669 bede 234 
‘Ttri’ 41204 14,545 abe 4137 # 435 abe 11,729 abc 30 ab 
‘Karim’ 3430 Pc 12,150 4 49.06 2 362 cd 8683 de az eho 
nN ‘Luiza’ 4129 @ 12,401 >ed = 42.35 bedef 346 4 9437 cde 30 ab 
‘Nassira’ 3659 abc 12,794 abcd 45.41 2 392 bed 8351 ¢ 24. cd 
Mean 3850 13,205 44.05 402 9574 27 
‘Faraj’ 3867 abe 14,900 # 41.48 Podet 464.2 12/104 3 2g abe 
‘Teri’ 3956 2 14,529 abc 42,39 bedef 454 ab 12,461 4 314 
as ‘Karim’ 3388 © 11,863 4 44,82 bed 359 cd 8713 de 29 ab 
‘Luiza’ 3843 abc 19235.4"6 39.08 f 400 abet 11,084 abcd 30a 
spe ‘Nassira’ 3545 be 14,777 ab 42.84 bedet 462 ab 11,156 abcd 25 bed 
Mean 3718 13,861 42.12 427 11,102 29 
‘Faraj’ 3564 be 13,274 abed 40.69 &f 448 ab 10,187 abcde a7 ed 
‘Teri’ 3904 abc 12,387 bed 39.13 f 402 abcd 10,992 abcd 314 
G3 ‘Karim’ 3689 abc 13,193 abcd 44,98 be £19 ehot Tipp set 31% 
‘Luiza’ 3851 abe 13,135 abcd 40.9 def 432 abe 11,084 abcd 304 
‘Nassira’ 3547 be 12,892 abcd = 44.49 bede 451 ab 9885 abcde 244 
Mean 3711 12,976 42.04 429 10,649 28 
Analysis of Variance 
Nitrogen Rate ns ns ns ns i ns 
Variety oe ne we + * we 
N Rate x Variety ns ns ns ns ns ns 


TKW = Thousand-kernel weight, Spk/m* = Spikes per square meter, G/m? = Grains per square meter, 
G/S = Grains per spike. ANOVA—Analysis of variance, where * = Significant at p < 0.05; ** = Significant 
at p < 0.01; *** = Significant at p < 0.001; ns = Not significant. Means without a common letter also vary greatly. 


The 2020 dataset provides critical insights into crop behavior under varying nitrogen 
conditions. Significantly, the ‘Faraj’ variety produced the largest yield, exceeding the overall 
mean yield of 1720 kg/ha. Its performance was particularly remarkable at the N2 level, 
where it registered a yield of 2178 kg/ha, surging ahead by 26.6% relative to the mean, as 
highlighted in Table 7. In stark contrast, ‘Karim’ and ‘Itri’ largely underperformed, with 
their yields languishing at approximately 22.0% and 14.7% below the mean, respectively. 
With respect to biomass accumulation, ‘Faraj’ produced the greatest biomass at the N1 and 
N2 levels, producing values, which were a substantial 33.5% and 27.2% higher than the 
average of 6305 kg/ha. ‘Karim’ was the least productive, with its biomass at the N3 level 
falling a concerning 25.6% below the mean. Intriguingly, while ‘Itri’ exhibited subpar yield 
and biomass figures, it was ‘Nassira’, which recorded the highest TKW of 42.41 g at the N3 
level—an increase of 8.5% above the mean TKW of 39.08 g. Concurrently, ‘Itri’ yielded less 
at the N2 level, with a TKW of 36.98 g—5.4% below the mean. 
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Table 5. Comparative study of wheat varieties under various nitrogen levels in 2018: yield, biomass 
and grain characteristics with ANOVA insights. 


Nitrogen Yield Biomass 


ee Rate VAnENy. (kg/ha) (kg/ha) TKW (g) Spk/m? Gim? G/s 
‘Faraj’ 4032 ¢4 10,443 <4 43,99 be 2654 7428 b Soh 
‘Ttri’ 4706 abc 13,993 @ 45.14 bc 3114 9723 ab 32 abe 
ae ‘Karim’ 4400 abcd 11,913 bed 46.9 ab 2924 8683 > 30 abc 
‘Luiza’ 4794 ab 14,009 4 43,92 bc 3214 9569 ab 31 abc 
‘Nassira’ 41224 12,682 bed 47.42 abe 311° 8308 > 27° 
Mean 4411 12,608 45.47 300 8743 30 
‘Faraj’ 427] abed 12,846 2b 46.01 abc 302 ° 9593 ab 47° 
‘Ttri’ 4124 abed 11,341 bed 43.46 © 3112 9266 » 31 abc 
N2 ‘Karim’ 4604 abcd 12,200 abc 46.75 2b 2812 8440 > 29 abe 
2018 ‘Luiza’ 4837 @ 11,397 bed 44.83 bc 3004 9520 ab 32 abc 
‘Nassira’ 4073 bed 10,827 bed 45.15 bc 284 4 6818 > 24¢ 
Mean 4382 11,722 45.24 294 8726 33 
‘Faraj’ 4293 abed 10,128 4 48.35 @ 2894 8238 > 29 be 
‘Teri’ 4576 abed 10,485 °4 44.78 be 3024 8940 > 30 abc 
NB ‘Karim’ 4763 abc 10,483 <4 46.58 abc 308 @ 12,679 @ 44 ab 
‘Luiza’ 4750 abc 11,186 bed 46.81 2 300 @ 9285 > 32 abe 
‘Nassira’ 4204 abed 10,422 <4 46.89 ab 2702 7441 > 30 abc 
Mean 4517 10,541 46.68 294 9315 33 
Analysis of Variance 
Nitrogen Rate ns cial * ns ns ns 
Variety ae ns ns ns i ns 
N Rate x Variety ns * ns ns ns ns 
TKW = Thousand-kernel weight, Spk/m* = Spikes per square meter, G/m” = Grains per square meter, 
G/S = Grains per spike. ANOVA—Analysis of variance, where * = Significant at p < 0.05; ** = Significant 
at p < 0.01; *** = Significant at p < 0.001; ns = Not significant. Means without a common letter also vary greatly. 
Table 6. Comparative study of wheat varieties under various nitrogen levels in 2019: yield, biomass 
and grain characteristics with ANOVA insights. 
Year oe Variety (aah) ‘hake TKW (g) Spk/m? G/m? G/S 
‘Faraj’ 3347 ab 11,165 2b¢ 39.69 ab 284 bed 7123 bed 22 abed 
‘Teri’ 2852 be 11,317 abc 36.54 4 254 bed 5989 <A opiebed 
N1 ‘Karim’ 2868 b¢ 9468 <4 39.42 ab 238 cd 5956 cd 22 bed 
‘Luiza’ 2853 be 10,287 bed 37.67 bed 246 bed 6070 <4 poibcd 
‘Nassira’ 2535 ab 8960 bed 39.74 230 bed 3864 bed 17 bed 
Mean 2891 10,239 38.6 251 5800 21 
2019 ‘Faraj’ 3863 4 12,475 ab 37.57 bed 367 4 12,1772 28 
‘Teri’ 2941 be 13,406? 35.36 de 302 abc 7790 bed 22 abed 
N2 ‘Karim’ 3262 > 10,635 Pe 38.04 bc 294 abed 8262 bc 24 abc 
‘Luiza’ 3294 ab 12,902 ab 37.31 bed 319 ab 9558 ab 25 ab 
‘Nassira’ 2879 be 10,225 <4 38.23 bc 267 bed 6488 bed 20 bed 


Mean 3248 11,929 37.3 311 8856 24 
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Table 6. Cont. 


Yield Biomass 


Year Rate Variety (kg/ha) (kg/ha) TKW (g) Spk/m? G/m2 G/S 
‘Faraj’ 3064 bc 9656 <4 39.61 a 275 bed 7190 bed 23 abed 
‘Ttri’ 2407 ¢ 9055 <4 33.82 ° 246 bed 5735 cd 19 <4 
N3 ‘Karim’ 2654. be 8742 <4 38.29 be 238 <4 6607 bcd 23 abcd 
‘Luiza’ 2.97 be 11,001 2b¢ 37.67 bed 267 bed 6561 bed 22 bed 
a0i8 ‘Nassira’ 2470 ¢ 7820 4 38.15 be 2274 4679 4 184 
Mean 2678 9255 37.51 251 6153 21 
Analysis of Variance 
Nitrogen Rate ia sacal al aad kang ns 
Variety - * Sap ns 7 ns 
N Rate x Variety ns ns ns ns ns ns 
TKW = Thousand-kernel weight, Spk/m* = Spikes per square meter, G/m? = Grains per square meter, 
G/S = Grains per spike. ANOVA—Analysis of variance, where * = Significant at p < 0.05; ** = Significant 
at p < 0.01; *** = Significant at p < 0.001; ns = Not significant. Means without a common letter also vary greatly. 
In the 2021 assessment—the year completing our longitudinal examination—the ‘Ttri’ 
variety displayed a complex but instructive performance profile. With a yield of 4319 kg/ha 
at the N1 nitrogen level, ‘Itri’ outperformed the overall mean yield of 3789 kg/ha by approx- 
imately 14.0%. However, it manifested a marked decline at the N3 nitrogen level, falling 
short of the mean by 11%. This dual behavior was also reflected in biomass accumulation; 
‘Itri’ recorded values 12.2% above the mean of 17,127 kg/ha at the N1 level but plummeted 
to 21.3% below the mean at the N3 level. Intriguingly, ‘Karim’ led in TKW at the N1 level, 
registering a TKW of 42.76 g, which was 6.5% above the overall mean of 40.14 g. In contrast, 
‘Itri’ marked the lowest TKW at both the N2 and N3 nitrogen levels, with values 7.4% and 
8.4% below the mean, respectively, as shown in Table 8. 
Table 7. Comparative study of wheat varieties under various nitrogen levels in 2020: yield, biomass 
and grain characteristics with ANOVA insights. 
Year ae Variety _Yield (kg/ha) ‘gle TKW (g) Spk/m2 Ghm? GIS 
‘Faraj’ 2079 2 84472 40.82 ab 208° 3777 abe 15 be 
‘Teri’ 1436 cde 5603 cde 38.06 cdef 135 bed 1968 4 13 Pc 
‘Karim’ 1323 ¢ 4733 ¢ 39.73 Pede 1164 19554 13° 
= ‘Luiza’ 1539 cde 6143 bede 37.36 ef 143 bed 2792 bed 14 be 
‘Nassira’ 2130 abc 7132 abcd 40.2 abc 194 abc 2930 bed 15 be 
Mean 1701 6411 39.23 159 2684 14 
‘Fara’ 2178 4 029 ab 40.29 abc 203 4 4069 a> 19a 
‘Ttri’ 1504 cde 5455 de 36.98 f 143 bed 2273 4 15 Pe 
N2 ‘Karim’ 1447 cde 5494 de 39.11 bedef 138 bed 22414 13° 
2020 ‘Luiza’ 1679 bede 6357 bede 37.08 f 162 abcd 2994 abed 15 Pe 
‘Nassira’ 1789 abc 6747 abed 39.82 bed 178 ab 2697 bed 13° 
Mean 1719 6416 38.66 165 2857 15 
‘Faraj’ 2022 a> 7581 abc 40.33 abc 208 4 4485 2 17 
‘Teri’ 1861 abc 6097 bede 37.15! 173 abc 3440 abed 17% 
N3 ‘Karim’ 1379 de 4692 ° 39.37 bedef 132 4 2408 <4 15 be 
‘Luiza’ 1685 bede 6015 cde 37.42 def 165 abc 2903 bed 15 abe 
‘Nassira’ 1754 abcd 6059 bede 42.414 165 abc 2819 bed 14 be 


Mean 1740 6089 39.34 167 3210 15 
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Table 7. Cont. 


Wass Nitrogen Variety Yield (kg/ha) aay TKW (g) Spk/m?2 Gim? G/S 
Analysis of Variance 
56 eee Rate ns ns ns ns ns ns 
Variety ro ro ce ee ro * 
N Rate x Variety ns ns ns ns ns ns 
TKW = Thousand-kernel weight, Spk/m* = Spikes per square meter, G/m? = Grains per square meter, 
G/S = Grains per spike. ANOVA—Analysis of variance, where * = Significant at p < 0.05; *** = Significant 
at p < 0.001; ns = Not significant. Means without a common letter also vary greatly. 
Table 8. Comparative study of wheat varieties under various nitrogen levels in 2021: yield, biomass 
and grain characteristics with ANOVA insights. 
meat ee vee (aah) ‘ko TKW (g) Spk/m? Gim? G/S 
‘Faraj’ 3887 abcd 18,989 @ 40.26 abed 413 abc 15,620 ab 30> 
‘Ttri’ 4319 @ 19,225 4 40.08 <4 405 abed 14,993 2 35 a 
Ni ‘Karim’ 4204 19,080 @ 42.764 421 ab 18,236 2 36 2 
‘Luiza’ 4203 ab 19,079 @ 39.99 cd 424? 19,227 39? 
‘Nassira’ 3544 abcd 14,634 abc 46.57 373 abcd 9204 ab 24a 
Mean 4031 18,201 41.93 408 15,455 33 
‘Faraj’ 3694 abed 15,857 abc 38.83 4 351 bed 13,030 2 34 ab 
‘Ttri’ 4124 abc 17,425 ab 37.18 © 378 abcd 14,685 2 38 ab 
N2 ‘Karim’ B78 eet 18,296 2 40.54 abed 410 abed 16,737 ab 32% 
‘Luiza’ 3556 bed 17,138 38.94 de B66 eee 16,419 2 33 2b 
2021 ‘Nassira’ 3956 abed 18,566 @ 41.88 abc 424° 18,382 ab 34% 
Mean 3813 17,456 39.48 392 15,852 34 
‘Faraj’ 3384 4 15,016 b¢ 39.21 de 348 <4 12,539 > 36 a> 
‘Ttri’ 3373 4 13,479 ¢ 36.75 © 3404 13,033 ab 35 ab 
N3 ‘Karim’ 3432 4 16,061 ®¢ 40.48 abcd 386 abed 16,745 > 65°» 
‘Luiza’ 3514 <4 14,946 Pe 38.42 de 343 cd 13,303 > 36 2 
‘Nassira’ 3915 abed 19,125 @ 40.23 bed 408 abcd 18,241 ab 37% 
Mean 3523 15,725 39.02 365 14,772 36 
Analysis of Variance 
Nitrogen Rate aes = pai id ns ns 
Variety ns ns ba ns ns ns 
N Rate x Variety ns ns ns ns ns ns 


TKW = Thousand-kernel weight, Spk/m* = Spikes per square meter, G/m? = Grains per square meter, 
G/S = Grains per spike. ANOVA—Analysis of variance, where ** = Significant at p < 0.01; *** = Significant 
at p < 0.001; ns = Not significant. Means without a common letter also vary greatly. 


The ANOVA results from the five-year period highlight the significant influence of 
wheat variety on the yield, biomass, TKW, Spk/ m2 and G/m?. The study also underscores 
the significant influence of nitrogen level on the yield, biomass, TKW and Spk/ m2 during 
the 2020/2021 growing season. This suggests that nitrogen level can significantly affect 
these traits, especially during growing seasons with greater rainfall, such as 2020/2021 
(429.8 mm). 
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3.5. Correlation Analysis: Yield, Yield Components and the Interaction between Nitrogen and 
Variety 

In this scientific investigation, we explored the intricate interconnections between crop 
yield and a multitude of agronomic variables, encompassing biomass, TKW, spikes per 
square meter (Spk/m7?) and grams per square meter (G/m?). The interrelationships among 
these variables are depicted in Figure 7, where the correlation coefficients, magnitude and 
orientation of each association are exhibited. 


Pearson's : 
Correlation Yield (q/ha) 
-10 -05 00 05 1.0 
= PM2 
-0.10 0.21 0.03 GIS 
ns 2 ns 
0.31 0.27 0.19 cM2 
om = ns 
as ee Bes on -0.03 = TKW@) 
ns ns ns ns ns 
> > 
C £ 
=. RS 
= 


ns P >= 0.05; * 7 < 0.05; ** ?< 0.01; and *** 7? < 0.001 


Figure 7. Pearson’s correlation matrix of key agronomic traits in wheat varieties under different 
nitrogen application rates. Note: This figure illustrates a Pearson’s correlation matrix for essential 
agronomic traits, namely the thousand-kernel weight (TKW), number of spikes per square meter 
(Spk/m7*), grains per square meter (G/m?) and number of grains per spike (GS). The significance of 
the correlations is indicated as follows: ‘ns’ denotes not significant; * denotes significance at p < 0.05; 
** denotes significance at p < 0.01; *** denotes significance at p < 0.001. 


In our comprehensive five-year study, we observed a strong correlation between crop 
yield and biomass (r = 0.60, p < 0.001), particularly in recently released wheat varieties, 
such as ‘Faraj’, ‘Luiza’ and ‘Itri’, which consistently outperformed older varieties in these 
aspects. Interestingly, while nitrogen application significantly enhanced yield and biomass, 
its impact on TKW and spikes per square meter (Spk/m7?) was less pronounced, suggesting 
a nuanced influence on different yield components. Furthermore, the uniform response of 
various genotypes to nitrogen fertilization indicates a stable interaction across varieties, 
underscoring the importance of genotype selection in tandem with tailored nitrogen man- 
agement for optimal wheat production. These findings highlight the complex interplay of 
genetic factors and agronomic practices in determining crop productivity, offering valuable 
insights for future agricultural strategies. 
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4. Discussion 


The findings of this investigation complement existing research, underscoring the 
pivotal role of both genotype selection and nitrogen application in influencing crop pro- 
ductivity. The marked superiority of recently released varieties, such as ‘Faraj’, ‘Luiza’ and 
‘Itri’, in terms of yield, biomass and yield components highlights the potential benefits of 
adopting these improved genotypes to boost crop productivity [25,26]. 

It was also observed that nitrogen application significantly influences the yield, 
biomass and TKW—a result, which aligns with prior studies reporting enhanced wheat pro- 
ductivity linked with optimal nitrogen fertilization [27-30]. Nevertheless, Spk/ m2, G/m?2 
and G/S largely demonstrated a non-response to nitrogen application; it is important to 
consider that the ANOVA analysis did reveal a near significant effect of nitrogen on Spk/m? 
(p = 0.071). This might hint at a developing trend not yet reaching conventional signifi- 
cance, suggesting that factors such as environmental conditions, soil fertility and specific 
management practices may exert a more substantial influence on these parameters [31]. 

One of the critical discoveries in this study was the lack of a significant interaction 
between nitrogen application and variety across all the investigated parameters. This out- 
come implies that the response of different genotypes to nitrogen fertilization is consistent 
across varieties. This finding is corroborated by previous research, such as that conducted 
by Peake et al. [32]. Their study revealed a noteworthy correlation among the genotype, 
environment and agronomic management with grain yield and lodging. This indicates 
that specific agronomic techniques, including particular nitrogen fertilization practices, 
could be broadly applicable across different varieties. Our study further emphasizes the 
importance of understanding the interplay of these variables in order to optimize grain 
yield. Moreover, an extensive examination by Kartseva et al. [33] revealed that diverse 
approaches to mitigating osmotic stress can be ascribed to variations in genetic make-up. 
This emphasizes that while the response to nitrogen fertilization may be generally uniform, 
the reaction to stress mitigation techniques can significantly vary across different genotypes 
due to their unique genetic characteristics [34]. 

Through a meticulous analysis of the association between crop yield and biomass, 
a robust and affirmative correlation was unearthed, thereby fortifying the prevailing 
understanding that the production of biomass plays a pivotal role in determining the 
ultimate grain yield. This aligns with the study by Balcha and Ayana [35], which found 
a positive correlation between grain yield and biomass yield in wheat genotypes. The 
importance of biomass in realizing high yield potential in wheat is a commonly echoed 
sentiment in earlier studies [36,37]. The observed positive correlations between yield and 
Spk/m?, G/m? and G/S indicate that enhancing plant density, grain production and the 
number of grains per spike could potentially lead to increased wheat yields [38-41]. 

An intriguing finding in this study is the negative correlation between TKW and 
Spk/ m2, suggesting a potential trade-off between kernel weight and plant density. This 
observation is corroborated by previous research [42,43]. The evidence of a trade-off 
between kernel weight and plant density further complicates the picture. Increasing 
plant density could potentially lead to smaller kernels, which could limit the potential 
yield increase from greater plant densities [44]. Hence, the development of varieties, 
which can maintain kernel weight while increasing plant density, may represent a fruitful 
area for future research and breeding efforts. The positive correlation between Spk/m? 
and G/m? indicates that an increase in plant density could potentially stimulate greater 
grain production. In contrast, the negative correlation between Spk/m? and G/S implies 
a potential reduction in the number of grains per spike as plant density increases. 

Despite the study revealing no significant effect of the interaction between nitrogen 
and variety on several parameters, it offers crucial insights. The consistent reaction of 
various wheat varieties to nitrogen fertilization among the genotypes is consistent with 
prior studies, which have shown comparable responses in wheat genotypes [45,46]. This 
consistency suggests that the nitrogen effect is independent of variety type. 


Nitrogen 2024, 5 


282 


However, the significant influence of nitrogen and variety as separate factors un- 
derscores the necessity for informed selection of the appropriate variety and nitrogen 
fertilization strategy to achieve large crop yields. Hence, this investigation underscores 
the significance of comprehensively comprehending the interplay between genetic and 
environmental factors in relation to wheat productivity, as this is crucial for augmenting 
crop productivity. 

Within the expansive framework of sustainable agriculture, the discoveries put forth in 
this study provide a persuasive argument in favor of adopting an integrated methodology 
for crop management. The successful cultivation of wheat requires a meticulous approach, 
which encompasses not only a careful choice of ideal wheat genotypes but also a profound 
comprehension of the intricate dynamics between various factors, including nitrogen 
fertilization, environmental conditions and management practices. With ever-increasing 
pressures on global food security due to population growth and climate change, such 
informed and targeted crop management strategies will be indispensable in achieving 
sustainable increases in crop productivity. 

Consequently, the results support the notion that breeders should strive to develop 
new varieties with improved nitrogen efficiency. A variety possessing enhanced nitrogen 
utilization capabilities has the potential to yield higher crop production while concurrently 
minimizing the requirement for fertilization, thus serving as a means to alleviate the eco- 
logical consequences associated with excessive nitrogen fertilizer application [47]. In light 
of the potential consequences associated with excessive application of nitrogen in agricul- 
tural systems, the aforementioned discoveries underscore the significance of conducting 
scientific investigations in order to foster the establishment of sustainable and ecologically 
sound farming methodologies. 

Another crucial aspect in this discussion is the relationship between plant density 
(Spk/m7?), grain production (G/m7*) and the number of grains per spike (G/S). The study’s 
findings suggest that a more thorough understanding of these relationships could offer 
the pathways for enhancing wheat yield. By identifying and implementing optimal plant 
densities, we may be able to boost grain production and, in turn, overall yield. However, 
caution should be exercised in pursuing this strategy, as there appears to be a negative 
trade-off between plant density and the number of grains per spike [48]. This observation 
is corroborated by previous research [49]. The evidence of a trade-off between kernel 
weight and plant density further complicates the picture. Increasing plant density could 
potentially lead to smaller kernels, which could limit the potential yield increase from 
greater plant densities. Hence, the development of varieties, which can maintain kernel 
weight while increasing plant density, may represent a fruitful area for future research and 
breeding efforts. 

In addition to the aforementioned observations, our comprehensive five-year inves- 
tigation (2017-2021) consistently demonstrated that diverse wheat varieties responded 
differently to varying nitrogen levels. Our study made the crucial observation that the 
annual climatic conditions were the main factor influencing these divergent responses. 
Temperature, rainfall, sunlight and other factors significantly affected the different varieties’ 
responses to nitrogen fertilization, suggesting a significant interaction between climatic 
conditions and nitrogen utilization efficiency. This underscores the need for adaptive 
management strategies, which account not only for the genetic attributes associated with 
diverse wheat varieties but also for the changing characteristics of the ecosystem in which 
they are grown, particularly in the newer varieties ‘Faraj’, ‘Luiza’ and ‘Itri’. Yield, biomass, 
TKW, the number of spikes per square meter (Spk/ m2), grains per square meter (G/ m7?) 
and the number of grains per spike (G/S) varied across the different wheat varieties and 
nitrogen rates, supporting the findings of Yadav et al. [50], who reported significant genetic 
variability in these traits among different wheat varieties. 

The ANOVA results over the five years indicate that both the variety and nitrogen 
rate significantly influence yield and other agronomic traits. This aligns with the findings 
of Sedri et al. [51], who reported that different physiological traits, including yield, were 
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significantly affected by various water conditions—an environmental factor comparable to 
nitrogen rate in our study. The precipitation data for each year also enrich the context of 
the results. For instance, 2018 saw the greatest precipitation and showed significant effects 
of nitrogen rate on biomass and TKW, possibly due to the increased water availability 
enhancing nitrogen’s impact on these traits. This is supported by Shibeshi and Kassa [52], 
who found that environmental conditions, including water availability, can significantly 
influence wheat yield and other agronomic traits. 

Additionally, the results suggest that the new varieties, ‘Luiza’ and ‘Itri’, have been 
optimized for high grain yield—a finding, which resonates with the complex trait interplay 
reported by Zhang et al. (2020) [53]. The ‘Faraj’ variety seems to have been bred with 
a focus on spike density—a trait highlighted for its importance in wheat yield by Joudi 
et al. (2017) [54]. The association of the ‘Nassira’ variety with high TKW suggests an 
emphasis on breeding for grain weight—a trait positively correlated with grain yield, as 
reported by Alonso et al. (2018) [55] and Racz et al. (2019) [56]. The positioning of the older 
variety, ‘Karim’, on the opposite side of the bi-plot suggests that it might possess different 
agronomic characteristics, likely due to different selection pressures during its development. 
These findings underscore the importance of targeted breeding strategies to optimize key 
agronomic traits, thereby enhancing wheat yield under various environmental conditions. 

The observed positive correlation between TKW and application of the greatest nitro- 
gen application rate (N1 120 U) is in line with the research conducted by Derkx et al. [57]. 
Their study on wheat demonstrated that increased nitrogen application rates resulted in 
increased nitrogen uptake, leading to the formation of larger—albeit fewer—grains. This 
phenomenon could be attributed to the presence of nitrogen, a vital element for plant 
physiological processes, which plays a significant role in increasing grain size and weight. 
The observations put forward by W6jtowicz et al. [58] provide additional evidence that the 
number of spikes per square meter significantly influences grain yield, especially when 
nitrogen doses are varied. The correlation between the density of spikes per unit area and 
G/m*, in conjunction with application of the N2 60 kg nitrogen application rate, suggests 
that a moderate level of nitrogen fertilization could potentially enhance these specific yield 
attributes. The findings by Guo et al. [59] support our current understanding of agronomic 
principles. They reported a noteworthy association between grain yield and two key factors: 
kernel number and thousand-kernel weight. This implies that specific clusters of wheat 
varieties have the potential to optimize yield traits, such as the number of grains per square 
meter. Moreover, the positioning of yield and biomass between the N1 and N2 nitrogen 
doses on the PCA plot is a significant finding. It suggests that an intermediate nitrogen 
dose, between N1 and N2, could potentially result in optimal yield and biomass. The 
findings of this study hold significant implications for the field of agronomy, particularly 
in relation to agricultural practices. The results indicate that a meticulous regulation of 
nitrogen application has the potential to enhance crop productivity to a considerable extent. 
This assertion is corroborated by the research conducted by Farkas et al. [60], wherein they 
discovered that wheat genotypes exhibiting favorable adaptability to diverse nitrogen lev- 
els exhibited the highest water use efficiency. Nevertheless, it is imperative to acknowledge 
that the aforementioned discoveries are predicated upon a single investigation, thereby 
necessitating additional scientific inquiry to validate and ascertain the precise nitrogen 
dosage, which optimizes both crop productivity and biomass accumulation. Future investi- 
gations should additionally take into account other variables, which have the potential to 
exert an influence on both crop yield and biomass production. These variables encompass, 
but are not limited to, soil composition, meteorological conditions and the specific genetic 
make-up of the cultivated crop. Huang et al. [61] suggested that good quality, high yield 
and high efficiency could be achieved through selection of high-quality wheat varieties and 
nitrogen accumulation and translocation. 

In conclusion, this comprehensive study underscores the importance of both the 
selection of an appropriate variety and optimal nitrogen fertilization strategy for wheat 
cultivation. However, the results also suggest that other variables, including plant pop- 
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ulation density and grain output, exert a substantial impact on crop yield and must be 
considered within a holistic agricultural management approach. Furthermore, the observed 
trade-offs among certain variables indicate that achieving optimal yield may necessitate a 
careful balance between these conflicting requirements. The collaboration between breed- 
ers and agronomists is crucial to optimizing the utilization of genetic advancements in 
conjunction with optimal agricultural techniques, thereby enhancing wheat productivity in 
a sustainable manner. In light of an escalating global food demand and the urgent need for 
environmentally friendly agricultural methods, the importance of adopting an integrated 
approach cannot be overstated. 

Production largely depends on the variety, year and their interaction, as indicated 
by Kolaric¢ et al., Popovic et al. and Stevanovi¢ et al. [62-64]. Profitable production based 
on high yields is impossible without the application of high doses of mineral or organic 
fertilizers, according to Stevanovi¢ et al. [64]. Fertilization is also important for grain qual- 
ity. These insights reinforce the complexity of agricultural systems and the multi-faceted 
approach needed for yield optimization. Studies of this nature yield invaluable insights, 
which can guide this undertaking, presenting pragmatic approaches to augmenting effi- 
ciency and resilience in the face of constantly mounting challenges. By cultivating a deeper 
understanding of the intricate dynamics between genetic predispositions and environmen- 
tal influences, we can enhance our preparedness to guarantee the sustenance of our food 
supply and advance the principles of sustainable agriculture in the forthcoming era. It is 
imperative to pursue further investigations in this domain, as augmenting our knowledge 
repository will enhance our capacity to effectively meet the escalating demands of our 
global population. Research endeavors such as this study serve as fundamental pillars 
in elucidating optimal methodologies and guiding subsequent investigations, thereby 
transcending the confines of wheat cultivation to encompass broader implications. The dis- 
coveries made by researchers significantly contribute to advancing our shared knowledge 
in the field of agronomy, with the ultimate goal of ensuring the long-term viability and 
improvement of global food security. 


5. Conclusions 


In our five-year investigation of Morocco’s Mediterranean climate, we systematically 
dissect the multi-faceted interactions between nitrogen fertilization, variety selection and 
annual variances, delineating their collective influence on wheat yield and the associated 
agronomic parameters. Our analyses underscore the significant yet nuanced impact of 
nitrogen application and varietal choice against the backdrop of year-to-year climatic 
fluctuations, revealing a pronounced variability in yield optimization strategies necessitated 
by these elements. The findings advocate for a more nuanced, adaptable approach to 
agronomic management, highlighting the necessity of integrating varietal selection and 
nitrogen application with an acute awareness of their interaction with seasonal climatic 
conditions. This approach is pivotal to enhancing wheat yield resilience and resource use 
efficiency, thus contributing to the broader objectives of sustainable agricultural practices. 
Our ongoing and future research will further elucidate the roles of plant density and 
comprehensive water management strategies, offering integrated insights for advancing 
agronomic efficiency and sustainability in wheat production. 
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